Role of Charged Gauge Fields in Generating Magnetic Seed Fields in Bubble 
Collisions during the Cosmological Electroweak Phase Transition 
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We calculate the magnetic field generated during bubble collisions in a first-order electroweak 
phase transition that may occur for some choices of parameters in the minimal supersymmetric 
Standard Model. We show that for sufficiently gentle collisions, where the Higgs field is relatively 
unperturbed in the bubble overlap region, the equations of motion can be linearized so that in 
the absence of fermions the charged W ± fields are the source of the electromagnetic current for 
generating the seed fields. Solutions of the equations of motion for the charged gauge fields and 
Maxwell's equations for the magnetic field in 0(1,2) space-time symmetry are expressed in closed 
form by applying boundary conditions at the time of collision. Our results indicate that the magnetic 
fields generated by charged fields in the collision are comparable to those found in previous work. 
The magnetic fields so produced could seed galactic and extra-galactic magnetic fields observed 
today. 
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I. INTRODUCTION 

Identifying the source of the observed large-scale 
galactic and extra-galactic magnetic fields remains 
an unresolved problem of astrophysics. One of the 
interesting possible sources is cosmological magneto- 
genesis, where the seed fields would have arisen dur- 
ing one of the early- universe phase transitions: the 
Quantum Chromodynamic chiral phase transition 
that led from the quark-gluon plasma to our present 
hadronic universe or the much earlier electroweak 
phase transition (EWPT) in which the Higgs and the 
other particles acquired their masses. See Ref. [1] for 
a comprehensive review. 

Here we consider magnetogenesis of the EWPT, 
which most likely requires a first-order phase tran- 
sition. Electroweak baroygenesis is a closely related 
problem of great interest and one requiring a first- 
order phase transition also. Although it is generally 



believed that there can be no first-order EWPT in 
the Standard Model there has been agreat deal 
of activity in supersymmetric extensions [3j , and for 
certain minimal extensions of the Standard Model 
there can be a first-order phase transition [1, d, 0] . 
One such extension that allows a first-order phase 
transition requires a right-handed Stop, the part- 
ner to the top quark, with a mass similar to the 
Higgs [1, 0|- Other models for CP violation and 
baryogenesis have been proposed, including two- 
Higgs models 0, @ and leptoquarks [|[ . 

Limits on parameter space of the minimal exten- 
sion of the Standard Model (MSSM) placed by elec- 
tric dipole moment measurements and dark matter 
searches allow a first-order electroweak phase transi- 
tion that could lead to successful electroweak baryo- 
genesis [9| and the possibility that we are explor- 
ing, namely that magnetic seed fields responsible 
for the large-scale magnetic fields seen today are 
created during the era of the EWPT. Accordingly, 
we consider the production mechanism of magnetic 
seed fields in a MSSM using an equation of mo- 
tion (EOM) approach similar to one we developed in 
studying bubble collisions in the QCD phase transi- 



2 



tion [nj. 

The most detailed previous research on the study 
of magnetic fields during the EWPT along these 
lines was based on an Abelian Higgs model [111 [l2l 
Il3j . in which the electromagnetic (em) field origi- 
nated from gradients in the Higgs phase when bub- 
bles of the broken phase collided and overlapped. In 
this paper, we calculate the magnetic field instead 
by starting from the MSSM using results of 0, [H| j 
where we derive the EOM directly from the elec- 
troweak (EW) Lagrangian and express the em cur- 
rent explicitly in terms of the charged fields ap- 
pearing in the Lagrangian. 

For cosmological magnetogenesis, the characteris- 
tics of the primordial magnetic field are vastly mod- 
ified during cosmic evolution to the present day. 
Much progress has been made in modeling this evo- 
lution in magnetohydronamics (MHD), first solving 
the equations for the period leading to the formation 
of galaxies and galactic clusters considering all rele- 
vant dissipative processes such as viscosity, and then 
for the evolution of these structures including the 
possibility that they provide a large-scale dynamo. 
Such studies have led to quantitative predictions for 
magnetic field energy and coherence length at the 
present epoch. One detailed study of possible galac- 
tic or extra-galactic magnetic fields evolving from 
the EWPT showed that random seed fields would 
not create such large-scale magn etic fields [16J. Re- 
sults of more recent studies [17(, however, are more 
optimistic and support the possibility that galactic 
cluster magnetic fields may in fact be entirely pri- 
mordial in origin. Our main interest in the present 
work is the mechanism that might have led to the 
creation of the seed fields during the EWPT rather 
than their cosmic evolution. 

For our present exploratory work, the precise na- 
ture of the extension of the Standard EW Model is 
not needed as long as it supports a first-order phase 
transition. We begin, in Sect. II, by reviewing the 
derivation of our EOM and explaining some of the 
considerations motivating our numerical work. In 
Sect. Ill we review the previous work done in the 
Abelian Higgs model. In our work as well as that 
of Ref. [ll|, [H, [Hj]j bubble nucleation and growth 
is driven by the dynamics of the Higgs field with 
an effective potential allowing a first-order phase 
transition, building on the classic work of Cole- 
man [l8l ]. Fermions are not explicitly considered. In 
Sect. IV we make use of the observations of Sect. II 
to simplify the theory along lines previously dis- 
cussed in Ref. 114], concentrating on a regime in 
which the bubble collisions may be considered " gen- 



tle". Restricting the application to gentle collisions 
has the advantage that the EOM linearize and dis- 
play a more transparent connection to the work of 
Elf E El- A closed-form solution of our EOM is 
given in Sect. V, and the results of our calculation of 
the corresponding seed fields are given in Sect. VI. 

II. EQUATION OF MOTION APPROACH 
TO MAGNETIC FIELD CREATION DURING 
THE EWPT 

Our ultimate goal is to determine the magnetic 
seed fields by solving EOM derived from the MSSM, 
a goal shared with Ref. [lj| . Both approaches obtain 
EOM derived from a similar Lagrangian, require a 
first-order phase transition, and take charged gauge 
fields as the source of the em current, but beyond 
this there are important distinctions. For example 
here we focus on a regime where linearized equations 
are appropriate, whereas in Ref. [la] the full nonlin- 
ear dependence on the W fields is retained assuming 
an I-spin ansatz for the W and for the em field. Ad- 
ditionally, here magnetic fields generated in bubble 
collisions are obtained, whereas in Ref. [l5j] em fields 
generated in nucleation are studied. Reference [la) 
explicitly includes the right-handed Stop in the La- 
grangian, whereas in the present work we integrate 
out all supersymmetric partners including the Stop. 

A. MSSM EW equations of motion 

In this section we summarize the derivation of our 
equations of motion [l4|] • As stated above, we imag- 
ine our MSSM Lagrangian to support a first-order 
phase transition, but all supersymmetric partners of 
the Standard Model fields arc integrated out, includ- 
ing the Stop. This Lagrangian is 

C MSSM = c l +£ 2 

+leptonic and quark interactions ( 1 ) 
£ 2 = \^- 9 -t-W,- 9 -B,)^-V(<S>), 

with 

W; v = dpWi-dvWi-gtijkWlWi (2) 
B^v = d^B v - d v B^ , 

where the W l , with i = (1,2), are the W ± fields, 
is the Higgs field, t % is the SU(2) generator. The 
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electromagnetic and Z fields are defined as 



a em 



1 



= (g'W 3 +gB^ 



= (gW 3 -g'B, 



(3) 



The effective Higgs potential V($) is not known at 
the present time, depending as it does on unknown 
parameters of the MSSM as well as the properties 
of the plasma in the early Universe at the time of 
the EW phase transition. The various parameters 
are discussed in many publications [|J. Fortunately, 
the specific form of V(&) is not relevant for the pur- 
poses of this paper, beyond the requirement that 
it should produce a first-order phase transition as 
it would in an underlying MSSM extension includ- 
ing a light right-handed Stop Q- Later, in an illus- 
trative calculation, we use a simple form for V{4>) 
taken from Ref. [lg]. For our calculations we need 
g = e/ sinOw — 0.646, g' — g tan#iy = 0.343, and 
G = gg'/^/j r T7 2 ~ = 0.303. 

In the picture we are using, the Higgs field plays 
a dynamic role in EW bubble nucleation and colli- 
sions. Our form for the Higgs field is 





Ax) 



(4) 



and thus 



We also use the definitions 

ftx) = p(x)e jQ(x) 
|0(x)| 2 = p{xf . 

With these definitions C 2 is (j = (1,2,3)) 

l\ 2 (wW J- <L\ 2 vt 2 



+ 3 '^)/2 + 0>[(|) 2 (^) 2 + { y -?B 
— ^-W 3 ■ B] — . 



(6) 



The equations of motion are obtained by minimiz- 
ing the action 



d*4C x +/: 2 ] = 0, 



(7) 



i.e., we do not include the leptonic parts of C. 
The result of doing this yields the following set of 



equations, where our metric is (1,-1,— 1,-1,— 1). 
The modulus p of the Higgs field satisfies the "p- 
equation" , 



= d 2 p(x) - i-p{x)[W 1 -W l 

+W 2 ■ W 2 } - p(x)ip v ip v + p(x) 



(8) 



dV 
dp 2 



where the quantity ip v is defined in terms of the 
phase of the Higgs field, the B field, and the W 3 
field as 



^ v ( x ) = d„e + g -B v - ^w 3 



(9) 



It obeys a "9-equation" , 

= d v p(x) 2 i>„(x), (10) 
with B field satisfying a "_B-equation" , 

= d 2 B„-d„d-B + p{x) 2 g'^(x). (11) 
Finally, the W field satisfies the two 'W-equations" , 

=d 2 W 3 -dvd-W 3 - f(x) 2 g^(x) (12) 

- ge 3 jk[W^d ■ W j + 2W j ■ dWj; - W^W k "} 
+ g 2 esjkWj i l e klm W l t l W™ 

for i — 3, and 

= d 2 W^-d u d-W i + ^f(x) 2 g 2 W l v (13) 

- ge ijk [W k d ■ W J + 2W J ■ dW k - W^W^] 
+ g 2 e iJ kW j > l 6 klm W lfi W™ 

for i = (1,2). 

These are exact equations of motion in our MSSM 
model. 



B. First-order electroweak phase transition 

As we indicated, in this paper we study the mag- 
netic field generated during a first-order electroweak 
phase transitions by finding solutions of the equa- 
tions of motion derived above in the particular 
regime of gentle collisions. Since the basis of our the- 
ory involves complicated coupled nonlinear partial 
differential equations (PDE), we would like to begin 
by illustrating some of the central ideas by discussing 
them within a model developed by Coleman [l8| that 
is conceptually transparent but contains much of the 
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complexity of the actual electroweak transition that 
we wish to examine. 

In Coleman's model, the Higgs field is a real scalar 
field <f> (9 = 0) satisfying the nonlinear PDE 



d 2 <t>{x) +cj)(x) 



dV 



= . 



(14) 



The effective potential V(4>), which specifies the de- 
pendence of the energy of vacuum on the scalar field, 
has a metastable minimum for a value of 4> — 4>i se P~ 
arated from a second minimum at (f> = (f>2 by a bar- 
rier. Coleman identifies the two minima with "true" 
and "false" states of the vacuum, respectively. In his 
model, a symmetry breaking term is added to V(<j)) 
to give the true vacuum a slightly lower energy. Ini- 
tially, the system is imagined to exist in the false 
vacuum for which ((f)) — <pi but over time makes a 
transition to the true vacuum in which (</>) = <p2- 

Coleman confirmed that in this model the system 
moves from the false to true vacuum by bubble nu- 
cleation as one expects of a first-order phase transi- 
tion. He showed that bubbles nucleate as tunneling 
(instanton) solutions of Eq. (fT4")) connecting (f> in the 
true and false vacua in Euclidean space-time, and he 
was able to calculate the nucleation probability. 

During the actual electroweak phase transition in 
the early Universe, the Higgs couples to the other 
particles in the electroweak theory, and as a result 
these particles acquire mass in the true vacuum as 
the value of (</>) changes from (j>x = to </>2 > 0. Once 
nucleated, bubbles grow and collide as solutions of 
Eq. (fT4"|) in Minkowski space-time. Coleman did not 
consider the physics of collisions, but these collisions 
of central interest in our work because it is through 
them that the magnetic field will be generated. 

We have modeled bubble collisions by obtaining 
numerical solutions of Eq. (| 14[) in 2 + 1 dimensions 
to illustrate an important feature of collisions. We 
assume that the bubbles nucleate at rest at time 
t = with radii r(t = 0) = Ro small compared 
to their radii r(t c ) = R c at the time of collision. 
For < t < t c the magnitude of the scalar field is 
constant at <j) = po = 5 (in arbitrary units) through- 
out each bubble, dropping rapidly to zero at its sur- 
face. Once the bubbles collide, the scalar field be- 
gins to fluctuate around the value po. The fluctua- 
tions are greatest for r ~ R c , just after the surfaces 
first touch, corresponding to the situation shown in 
Fig. [TJ A vertical slice through the symmetry axis 
is shown in Fig. O It is seen that at this point 
the scalar field fluctuates by as much as 30% from 
p = pa . These fluctuations rapidly damp out and es- 
sentially disappear once the bubbles interpenetrate 




FIG. 1: Two bubbles colliding in the Coleman model as 
explained in the text. The distance scale is in arbitrary 
units. 
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FIG. 2: A vertical slice through the colliding bubbles 
shown in Fig. 1 as explained in the text. The distance 
scale is in arbitrary units. 



with r(t) « R Q + R c . 

For the purposes of this paper it is important that 
in the collision these fluctuations in the scalar field 
remain small compared to po- We refer to a colli- 
sion of this character as a gentle collisions. We will 
make use of the observation that collisions between 
bubbles tend to be gentle to guide our approach to 
show how magnetic field generation arises as a result 
of charged dynamics as the bubbles collide. 

To model the dynamics of a first-order electroweak 
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phase transition of the early Universe, one needs to 
consider the dependence of the effective potential 
V($,T) on temperature T (see Ref. for a com- 
prehensive phenomenological study of the kinetics 
of cosmological first-order phase transitions such as 
the EW phase transition in terms of such an effec- 
tive potential). Just as imagined in Coleman's toy 
model, the vacuum state of the Universe corresponds 
to a local minimum in U(<£>, T). For T greater than a 
critical temperature T c , the minimum in the effective 
potential is at ($) = 0. At T c , V(<&,T) develops a 
degenerate second minimum at a value (<£>) > sep- 
arated from the firsr by a barrier. As the Universe 
continues to expand and cool, the second minimum 
becomes deeper, meaning that the Universe begins 
to tunnel from the original, now metastable, false 
vacuum to the true vacuum by bubble nucleation. 
In Ref. T c is estimated to be T c « 166 GeV. 
Because the two minima are separated by a barrier, 
the transition is delayed, a process referred to as su- 
percooling. As the phase transition develops, the 
bubbles expand with a speed v and eventually com- 
pletely merge, at which point the phase transition 
is completed. The terminal speed depends on the 
parameters of the phase transition but appears to 
be quite low, perhaps v wa u w .lc, as discussed in 
Ref. [H. 

Our theory is more complicated than that of Cole- 
man since it based on the MSSM, with the Higgs 
coupled to the other fields of that theory. We note 
here, by reference to Eq. |5]), that the gentle charac- 
ter of the collision illustrated in Fig. [5] may change 
when the W fields become large, with the Higgs field 
showing strong variation in the collision region. We 
refer to collisions of this character as violent colli- 
sions, the subject of a future publication. 



III. BUBBLE COLLISIONS IN THE 
ABELIAN HIGGS MODEL 



The Abelian Higgs model has been of interest as 
a prototype for the generation of magnetic fields in 
the early universe in collisions of bubbles during a 
first-order EW phase transition [ll|, Q3, Gj|. The 
Lagrangian of the Abelian Higgs model describes a 
complex scalar field coupled to the em field A^. It 
corresponds to the Lagrangian Eq. (JTJ) in the abelian 
sector, formed by eliminating the W fields, identify- 
ing A^ 71 with the field B^, and relating the electric 
charge e to coupling parameter g' as e = g' /2. In 
this limit, the equations of motion for the Higgs and 
magnetic field may be read off the results in Eqs. |T 



EE), 



and 



d v p(x) 2 Mx)=0 , (15) 

d 2 B v - d^d v B„ + 2ep(x) 2 ^(x) = 0, (16) 

8V 

d 2 p(x)-p(x)ip v ^+p{x)—=0. (17) 



dp 2 



The quantity tpv is now 



(18) 



In this case, the phase transition is driven by the 
dynamics of the Higgs field along the lines studied 
by Coleman [l8[ and discussed in connection with 
Fig. Q] and Fig. [2] in the previous section. 



A. Gentle bubble collisions in the Abelian 
Higgs model 

In their analysis of the Abelian Higgs model, Kib- 
ble and Vilcnkin [Til ] suggested one way in which 
magnetic fields might be generated as bubbles col- 
lide. They considered the regime of gentle collisions, 
where p(x) remains constant, or nearly constant, in 
the region of overlap of the colliding bubbles. They 
were able to gain insight into the generation of mag- 
netic fields in this case by making an expansion 
about point p(x) — po, which we shall refer to as 
the Kibble- Vilenkin point. For the case of gentle 
collisions, we can assume the following expansion, 



p(x) 



= pa + aSp(x) 

= ai/> (1) + a 2 ^ 2) + 



(19) 



where a is the magnitude of po — p(x) and is small 
by assumption (refer to the discussion in connection 
with Fig. [1] and Fig. ^ . Substituting the expansion 
into the equations of motion and requiring that the 
equations be satisfied at each order in the expansion 
parameter a, the relevant and B equations give, 
to leading order in a, the following results 



where now 



(d 2 + 2e 2 pM»=0 
d>« = , 

= d v Q + eAt m W 



(20) 
(21) 

(22) 



From these equations one can easily derive the equa- 
tions of Kibble and Vilenkin [ll| in the axial gauge 



for A" n . 



G 



B. Bubble collisions and 0(1, 2) symmetry 

Prior to nucleation, the dynamics is most easily 
formulated in the Euclidean metric, which is 0(4) 
symmetric. After nucleation, the bubble expands 
in Minkowski space-time, which is 0(1,3) symmet- 
ric For collisions, Kibble and Vilenkin study the 
case of two bubbles nucleated simultaneously at time 
t = at positions symmetrically located about the 
origin on the z— axis. In the notation of Ref. (l3| . 
the radii of the bubble at nucleation are Rq. Once 
nucleated, the bubbles expand with radii r(t) ac- 
cording to 

r{tf = R 2 + t 2 . (23) 

The bubbles collide on the z axis at t = t c with radii 
r(t c ) — R c , and hence the nucleation points are at 
z = ±R C . In this system, the collision of two bubbles 
has the symmetry 0(1,2) and may be expressed in 
the coordinates (r, z), where 

t 2 = t 2 - x 2 - y 2 . (24) 

Kibble and Vilenkin formulated their Abelian 
Higgs model in terms of the variables (r, z) and 
found the 0(1,2) symmetric solutions of Eqs. (|2H|) 
and ([22)) . They demonstrated from them that when 
the phase of the Higgs fields is initially different 
within each bubble an axial magnetic field forms as 
the bubbles collide and that this field has the struc- 
ture of an expanding ring encircling the overlap re- 
gion of the colliding bubbles. 

C. Finite wall speed and plasma conductivity 

The Abelian Higgs model has the feature that the 
bubbles accelerate freely to the speed of light along 
space-time hyperboloids. Corrections come from a 
variety of sources and were discussed comprehen- 
sively in Refs. [ll|, [l^, HH, an d the effects of these 
corrections on the subsequent evolution of the mag- 
netic field was estimated. We discuss these results 
briefly next. 

One of the important corrections is that bubble 
walls reach a terminal speed v wa u < < 1 as the bub- 
bles experience collisions with constituents of the 
plasma. Another important correction is the con- 
ductivity of the medium. The effects of finite con- 
ductivity lead to decay of the currents (and there- 
fore the ma gne tic field) with a characteristic time 
td ~ cr/m 2 [111 ]. However, with values of conduc- 
tivity that are believed to characterize the plasma, 



currents and magnetic fields persist on time scales 
that are long compared to those of the symmetry 
breaking scale. 

Another consequence of the large conductivity 
arises from the following considerations. Since the 
magnetic fields propagate with the speed of light, 
in the absence of conductivity, for slowly expand- 
ing bubbles these fields would very quickly escape 
from the region of the bubble collision and move into 
the surrounding false vacuum. However, because of 
the large conductivity the magnetic fields become 
"frozen" or confined to the interior of the bubbles. 
This is particularly important in the case of slowly 
moving bubble walls, where this effect prevents the 
escape of magnetic flux into the false vacuum. Kib- 
ble and Vilenkin showed that the loss of flux is neg- 
ligible provided that aR c v » 1, where R c is the 
bubble radius at collision time. 



IV. MAGNETIC FIELD GENERATION IN 
THE MSSM IN GENTLE COLLISIONS 

The main difference of our current work compared 
to the Abelian Higgs model arises from the fact that 
in our MSSM model the source of the current is the 
charged gauge fields. These IV ± originate in the 
plasma and populate the bubble by passing through 
the bubble wall following nucleation. As they do 
this, they acquire mass and thereby loose kinetic 
energy. We envision them to cool as the bubble ex- 
pands, thereby occupying a low-energy mode, the 
occupation of which grows as the bubble expands. 
We assume in this work that early in the phase tran- 
sition all such entering the bubble drop into this 
mode, which has a high degree of coherence, much 
like a state of electrons in a superconductor (except 
that the W are bosons) . It is quite different from the 
more familiar thermal modes, which are incoherent. 

As we shall see, this mode plays a special role in 
the theory. A significant consequence of the strong 
coupling to the Higgs field is that the W fields follow 
the evolution of the Higgs field, and this is described 
by our EOM. We solve our EOM the magnetic field 
in gentle bubble collisions for the MSSM in paral- 
lel to the Abelian Higgs model by applying simple 
boundary conditions in Kibble- Vilenkin geometry at 
the time of collision. From our results we will be able 
to ascertain the relative importance of the charged 
in the MSSM as compared to the Higgs mech- 
anism of the Abelian Higgs model as a source of 
magnetic field generation. 

The dynamics of the coupled fields is clearly a non- 
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equilibrium component of the phase transition, and 
eventually the coherence must dissipate as the sys- 
tem thermalizes. Determining the value time scale 
of the dissipation is an important issue that we do 
not address here. We assume that the coherent evo- 
lution of these fields after the collision give rise to 
a magnetic field before thermal equilibrium is re- 
established. 

It is our eventual goal to determine the magnetic 
field collisions between individual bubbles possess- 
ing 0(3) symmetry as they would when evolving 
from nucleation with a finite wall velocity. Thus, 
in the future we will examine the more numerically 
intensive case of collisions by solving the equations 
of motion numerically. Numerical methods for deal- 
ing this situation will be reported in forthcoming 
publications [2l| . 

Here, we first derive, in Sect. IV. A, the EOM 
in the case of gentle collisions. In the MSSM, the 
EOM are complicated nonlinear PDE coupling the 
W, B, and <fi (Higgs) fields. From the solution of 
these equations, the physical Z and A em fields are 
determined by Eq. (J3]). As shown in Sect. IV. B, this 
leads in turn to an expression for the electromagnetic 
current in terms of the and to a corresponding 
Maxwell equation. 



A. EOM for gentle collisions in electroweak 
theory 

As indicated, we consider here only the regime of 
gentle collisions in the context of the full MSSM. 
This leads to a linearized theory giving EOM for 
these fields. Although the linearization leads to the 
same limitations as the Abelian Higgs model dis- 
cussed at the end of the last section, fortunately the 
corrections for the finite wall propagation and the fi- 
nite conductivity of the plasma are easily estimated 
since the considerations are similar to those of the 
Abelian Higgs model. 

As in the Abelian Higgs model, the MSSM equa- 
tions simplify upon expansion about the Kibble- 
Vilenkin point, 

p{x) = po + aSp(x). (25) 

The fact that tp and W d (for d = (1,2)) enter 
quadratically in the p-equation places two impor- 
tant constraints on these quantities: (1) ip an d W d 
must have an expansion in odd powers of a 1 / 2 , if we 
require the square of these quantities be analytic in 
a; and, (2) expanding this equation to leading order 
in a 1 / 2 , we find that the terms V> , w^ 1 , and w^ 2 



must vanish. This is most easily seen in the Eu- 
clidean metric, from the fact that the square of each 
enters with the same sign. However, the same must 
be true in the Minkowski metric as well by analytic 
continuation. In view of these considerations, ip„ 
and W d for d = (1,2) have the following expansion 

i( I )=a 1 /V« + a 3 / 2 #+ ... (26) 

W d = a^ 2 w^ d + a^ 2 w^ d + ... . (27) 

It is natural that an expansion in the same param- 
eter a 1 / 2 remains appropriate for d = 3. However, 
there is no requirement that the leading term vanish, 
so we take 

W 3 = w^ 3 + a l ' 2 w^ + a 3 / 2 w^ 3 + ... . (28) 

The B-, 0-, and VF-equations then give, to first 
order in a 1 / 2 , 

[9 2 + f(5 2 +5 /2 M X) =0 (29) 

d>« - (30) 

where now 

^\x) = d a e- {g +9 2 }/ zw. (si) 

Comparing these equations to those in the Abelian 
Higgs model, Eqs. (|20I21I22|) . we see that the Z field 
here plays the same role as the em field did in that 
case. Specifically, in this case, the phase difference 
of the Higgs field now determines the gauge field Z^ 
of Eq.([3]) within the bubble overlap region. Thus, 
for gentle collisions, the mathematical problem in 
leading order is no more complicated that it was in 
the abelian case. 

Equations for wip d may be obtained by expanding 
the B- and VF-equations through order a 1 / 2 . For d — 
1 or 2 (corresponding to d' = 2 or 1, respectively), 
we obtain the pair of equations 

= d 2 w^ d - d v d ■ w (1)d + m 2 w^ d (32) 

- 2[d^(w^ 3 w^ d ' ~ w^ d 'w^ 3 ) 
+ (w^ d 'd^ 3 -w^ 3 d^ d ') 

- (w^ d 'd v w^ 3 - w^dM 1),xd ')] 

- 4[(w (Q)3 ) 2 w^ d - w^ 3 ■ w il)d w^ 3 } , 

where m is the mass of the W field and is given by 
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The corresponding equation determining w v for 
d = 3 is 

d 2 w^ 3 - d u d ■ w^ 3 = plg^ , (34) 

which can be solved once the driving term ip^'(x) 
has been independently determined from the solu- 
tion of Eqs. p9l30]) . 

Note that knowledge of the field wi for d = 3 
is required in order to solve Eq. (f32|) . This field is 
found to be the solution of 

3 2 w^ 3 ~ d v d ■ w^ 3 = . (35) 

Because the mass term is absent in this equation, 
clearly «/°) 3 can not contribute to Eq. (|28|) inside the 
bubble where m is non-zero. Thus w^ 3 must vanish 
inside the bubble. Consequently, the equations for 
-u^ 1 -* 1 and u^ 1 ) 2 simplify and become 

= d 2 w a v - d v d ■ w a + m 2 w a u , (36) 

and we see that for sufficiently gentle collisions, all 
relevant equations are linear inside the bubble. 



B. The electromagnetic current for gentle 
collisions 

We may find the Maxwell equation for the elec- 
tromagnetic field A e J n {x) by taking the linear com- 
bination of the W^ 1 and and B indicated in Eq ©. 
Using Eqs. ( [TTj) . (fT!?l ). an expression for the corre- 
sponding em current j„(x) immediately follows. It 
consists of terms quadratic and cubic in the three 
fields W l {x). 

This result for j v {x) may again be simplified by 
expanding the A em and W fields in powers of a 1 / 2 . 
Letting (x) refer to the terms in the expansion of 
A e J a {x), we find that the leading term of A e J n (x) is 
a v (x), satisfying the following Maxwell Equation, 

d 2 a^ ~ dvd -of® (37) 

- w^ a d v w^ b + 2w^ a ■ dw^ b ) 
= 4irj„(x) , 

From this we learn that the first non-vanishing con- 
tribution to the em current is of order a 3 / 2 and that 
it depends on the components Wv of the charged 
W l fields (i — 1 and 2), calculated at order a 1 / 2 . 



Expressing the current in terms of w l of Eq. (134j) . 
we find 



„„',ab3 

- w^ a d v w^ b + 2w^ a ■ dw^ b ) , 

It is easy to prove that this current is conserved, 

d-j(x)=Q, (39) 

since w^(x) and w b (x) that appear in Eq. p8[) com- 
mute with each other and satisfy the equation of 
motion given in Eq. (|36|) . 

One of the most important features of the deriva- 
tion is that the source current of the electromag- 
netic field is given by the charged gauge VF ± fields. 
This is expected physicall y, a nd is in sharp contrast 
with the Kibble- Vilenkin Till. Ahonen-Enqvist [T3] , 
Copeland-Saffin-Tornkvist [131 ] picture, in which the 
source for the electromagnetic field arises entirely 
from the Higgs. 

It should be clear from the fact that the electro- 
magnetic current in Eq. Q38p is antisymmetric in the 
labels a and b that this current will vanish unless 
the field W a (x) for a = 1 has a different dependence 
on x from that for a = 2. Thus, in an isolated bub- 
ble the electromagnetic current will vanish if W 1 (x) 
and W 2 (x) in this bubble differ at most by a phase. 
However, when two bubbles meeting this condition 
collide, it is in general not the case that W a (x) satis- 
fying the equations of motion will be proportional in 
the region of overlap. This is the reason why mag- 
netic fields arc in general produced when bubbles 
collide. We model this below by solving the equa- 
tions for W a with different sets of boundary condi- 
tions on W a for a — 1 and a = 2. 



V. SOLUTION OF THE MSSM EQUATIONS 
OF MOTION FOR BUBBLE COLLISIONS IN 
0(1,2) SYMMETRY 

Following Kibble and Vilenkin, in Sect. V.A we ex- 
press the EOM using the (r, z) variables and find the 
general solution of the EOM for fields for a pair 
of colliding bubbles. In Sect. V.B we find the solu- 
tions for specific boundary conditions. In Sect. V.C 
we determine the current corresponding to these so- 
lutions and find expressions for the corresponding 
magnetic field by solving Maxwell's equations. 
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A. W equations with 0(1,2) symmetry 



B. Boundary conditions 



To express the equations of motion, Eq. (|36|) , in 
terms of the (r, z) coordinates we define 

wl{x) = <(t,z), v = i 

w a v {x) = x v w a {T,z), v = (0,1,2) , (40) 

with w" = —w za . (In the remainder of the paper, we 
will find it convenient to use a to denote the Lorentz 
index for the values v = (0, 1, 2).) We then find 



(^ 
K dr 2 



9 r) r) 2 

^-^ + mV(x) = 

T OT OZ l 



and 



dr 2 



4 d 
t dr 



if 

dz : 



+ m 2 )w a {x) = 



where we have used the fact that 

dw za „ „ dw a 



d-w a (x) 



3w a 



= 



(41) 



(42) 



(43) 



dz dr 

Equation (f43|) is a consequence of Eq. (|36|) and fol- 
lows by taking the four-divergence of this equation. 
Comparing our results to the Abelian Higgs model, 
we see that w a takes over the role of the electromag- 
netic vector potential and d z w za takes over the role 
of —m/e<d. In our theory, to find the electromag- 
netic field we have the additional steps of construct- 
ing the electromagnetic current in Eq. (f3"5| and then 
solving Maxwell's equation, Eq. ([571) . 

The solution of Eqs. (|41I421 and O is found by 
standard methods. Expressing W"(x) as a Fourier 
transform in z, Eqs. (|41l42p give ordinary differential 
equations for the r-dependence, yielding 



w a {r,z) = i,/- 

a k (coswfcT- 
b k (sinWftT + 



= ( 



sinwfer, 

U! k T ' 
COSWfcT 

WfcT 



(44) 



))dk 



and 



if!" 



r dkz 



(cfeSinwfcT (45) 



-I- dfc cos io k r)dk 



We may use the auxiliary condition of Eq. (|43j) to re- 
late the coefficients of the linearly independent func- 
tions in the expansion of w a and w za . This gives 

ic k = a k uj k /k , (46) 
id k = b k uj k /k . (47) 



and 



As we indicated earlier, the relative phase of the 
two charged fields is irrelevant for a single bub- 
ble, but in a collision, this is not so. To be specific, 
let us assume that when the colliding bubbles first 
touch at z = 0, the magnitude of the fields is 
the same in each bubble. It will sometimes happen, 
though, that the sign of Wg is opposite in the two 
colliding bubbles while W~ the same sign. We will 
refer to the first situation, where W z changes sign 
across the point of contact at z = 0, with as super- 
script I, and the second situation with a superscript 
II. 

The two situations correspond to distinct solu- 
tions of the EOM, each described by its own bound- 
ary condition. The boundary condition BCI for the 
first situation is 



w zl (t = t c , z) = w e(z) 



(48) 



Ot 



w z1 {t = t c ,z) 







where e(z) is the sign of 2, and that for the second 
situation BCII, 



w ZI1 (r = t c ,z) 



(49) 



,.zll 



dr 



(r = t c ,z) = . 



We note in passing that, from the definition of r in 
Eq. ((24)) . the boundary conditions in Eqs. (|48]) and 
specify the value of w z on a cylinder of radius 



\J x 1 + y 2 = ^Jt 2 — t 2 = bit), where b(t) is the radius 
of the circle of intersection of the colliding bubbles 
at time t. 

The constant w in Eqs. (|48l and |49|) fixes the den- 
sity of the charged gauge fields W ± inside the bub- 
ble; charge neutrality requires that w is the same 
for both a — 1 and a — 2. We determine w in 
the appendix in a model that assumes the average 
number density of the W ± quanta inside the bubble 
(in the true vacuum) is constant in time and equal 
to the number density of those quanta in the 
thermal plasma that can make a transition into the 
bubble without violating overall energy-momentum 
conservation. Since the charged gauge field inside 
the bubble is normalized to the one outside the bub- 
ble, we should identify in Eq. (f2"Tj) with w^ d 
rather than with a 1 ^ 2 w^ d . 

From Eqs. (|44I45|) we can determine the full (t, z) 
dependence of the fields w za and w a corresponding 
to each of the two boundary conditions. We will then 
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calculate the magnetic field from a collision of the 
two bubbles as an indication of what magnetic fields 
can be expected when they arise from the current of 
the charged W ± fields. 

Equating Eq. (|45p and its r derivative at r = t c 
provides, in a straightforward fashion, the coeffi- 
cients a,t and bk for BCI and BCII. The solution 
is most direct for BCII, leading to 



w z (r, z) = (sin mT + mt c cos mT) (50) 

mT 

w h (t,z) = , 

where T = r — t c . For BCI, we may perform the sum 
over the modes k to obtain 



wt, 



w zI ( T ,z) = — e(z)e(T-\z\ 



(51) 



( y + A)j ( mv /T 2 - z l2 dz' 
+ 0(|z| _T)l(i + A) sin mT 
7 (r,z) = -^Q(T - \z\)(J Q (my/f^z^) 



yr 2 - z 2 

mt c T 



Jiimy/T 2 - z 2 )) 



and 



4nj z (r,z) = -—$ = ={2tw 1 ^- (54) 



VS 2 + .9' 2 

dw zI 



Or 



w 



zll 



dz 



When evaluating the partial derivatives in 
Eq. (|53I54|) using the expressions in Eqs. (fBTil 
l5"Tj) we do not let the derivatives act on the 9 func- 
tions. We note that ignoring the surface derivatives 
is quite consistent with current conservation and 
should therefore lead to a valid expression for the 
magnetic field throughout the bubble interior. 

Because the electromagnetic current has the form 
given in Eq. ([52]) . the electromagnetic field has this 
form also, 



a v {r,z) = (a z (T,z),x a a(T,z)) 



(55) 



Maxwell's equation becomes quite simple in the 
(r, z) with the axial gauge, a z = 0, namely, 



dz 



^a(r,z) = 4tt/(x,t) . 



(56) 



Applying the boundary conditions on the a(r, z) 
field, namely a(t c , z) = and d z a(r = 0, z), we find 



C. Magnetic field 



a(r,z) = -4w dz' j(T,z")dz" . (57) 



To obtain the magnetic field we need to solve 
Maxwell's equation, Eq. ([37]), with the current given 
by Eq. (|3"5)) . and with the fields appearing in 
the current given by the solutions of Eqs. (|4T1 14*21 
and 23]). The solution of Maxwell's equation is com- 
pletely determined once we specify the boundary 
conditions on the A em field with the fields given 
above in Eqs. ([SD I IST ]) . 

In the (t, z) coordinates, the current in Eq. 
may be written in the form 



The magnetic field B l is the curl of A 1 em where 
A 1 em — (xa(r, z),ya(T, z),0) in the axial gauge with 
a(r,z) given in Eq. ([BT]) . Thus, 



B z = 0, 



B x = -47TJ/ f j(r,z')dz' 
B v = 4ttx / j(r,z')dz' . 

J —OO 



(58) 



3u{x) = (jz(T,z),X a j(T, z)) 



(52) 



Evaluating the current with fields satisfying the 
boundary conditions given in Eqs. (H2] and 05]) and 
using the auxiliary condition of Eq. (143[) . we find 
that Eq. (|3"5|) simplifies and becomes 



Clearly, as in the work of Kibble and Vilenkin, the 
B field is entirely azimuthal, encircling the axis of 
collision. 



VI. NUMERICAL RESULTS 



47rj(r,z) = 



qq' 1 , ,jt dw 1 

-(-2tw zI1 — (53) 



zI dw 



12 T 



zll 



w 



zll 



Ot 



dw zI 
8t 



dz 
) 



In this section we compute the magnetic field us- 
ing the theory developed above, which assumes a 
non-conducting medium a = and a terminal wall 
speed v wa ii = 1- We will compare our magnetic 
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fields to the results obtained in the Abelian Higgs 
model under the same assumptions. 

Although, as we admitted earlier, the assumptions 
of (7 — and v wa ii = 1 are unrealistic for the ac- 
tual EWPT, the corrections have been thoroughly 
studied in Refs. [ll| and [12j |. The corrections are 
therefore well-understood, and they will be easy to 
apply to our results since they are not specific to the 
source of the electromagnetic current. We will also 
do this below. 

We first show our azimuthal field using our 
results in Eqs. (gEl 1491531 and [551 )• The mass ap- 
pearing in the expressions for the W fields is of 
course the mass of the W boson, m\y = 80.6 GeV. 
To facilitate the comparison to the Abelian Higgs 
model in Ref. [l^], we will assume, as there, that 
bubbles nucleate at points symmetrically placed on 
the z-axis at ±i? at time t = 0, that they expand 
from the point of nucleation with a speed approach- 
ing the speed of light, and that they collide at time 
t = t c = R. They present their results in terms of 
the time after collision, which we will call St, so that 
t = tc + St and r 



^(St + R) 2 - p 2 . 
We write the overall normalization of the electro- 
magnetic current by combining the factor contain- 
ing the coupling parameters appearing explicitly in 
Eq. (|53[) and the square of the normalization of the 

fields appearing explicitly in Eqs. (1481 and [ 
This normalization is thus 



w 



99 



Vi 2 



(59) 



Using the values of g and g' quoted below Eq. ^ and 
the value of w specified in Eq. (f6"5)) of the Appendix, 
we find that Eq. (|5"9"]) becomes 



w 



99' 



V9 2 + f 



w -38.5+ 1.36T C GeV (60) 
= 2.32mw , 



taking T c = 166GeV from Ref. [l9(. 

Figure. [3] shows our calculated field at a sequence 
of times St after the time of collision assuming that 
the bubbles nucleated at a distance of ±10 from the 
origin (R = 10). The magnetic field shown in this 
figure is the value of the azimuthal field in the sym- 
metry plane, which is the plane perpendicular to the 
axis of collision at the point where the bubbles first 
collide. Distance and time are expressed in units of 
l/m w and the magnetic field in units of m^. It 
is clear that the strength of the field is largest in 
the outermost region of the surface of the expanding 
overlap region and that this region narrows as the 
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FIG. 3: Magnitude of the azimuthal magnetic field cal- 
culated in the theory of the present paper. The field is 
shown as a function of distance p from the axis of colli- 
sion in the symmetry plane at time St — 5, 10, 15, and 
20. 
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FIG. 4: Magnitude of the azimuthal magnetic field cal- 
culated in the Abelian Higgs model. The field is shown 
as a function of distance p from the axis of collision in 
the symmetry plane at time St — 5, 10, 15, and 20. 



bubble expands. The fact that the speed of expan- 
sion can be superluminal is not in contradiction with 
relativity, as discussed in Ref. [l3| . 

This is to be compared to the result in the Abelian 
Higgs model, shown in Fig. [4] A convenient expres- 
sion for the magnetic field is given in Eq. (29) of 
Ref. pj]. We have taken [l| 6> = 1, the charge 
e = 1, again expressing distance and time in units of 
inverse mass and the magnetic field in mass squared, 
where the mass in this case is that of the vector bo- 
son in that model. 

There are several points to be made from the com- 
parison of Fig. [3] and Fig. [4] The first point is that 
the magnitude of the magnetic field in the present 
model is about the same size as the one obtained 
in the Abelian Higgs model. This is the case even 
though the description of the physics is quite dif- 
ferent in the two cases, with the source of the em 
current in the present model being charged gauge 
bosons that originate in the plasma. Once these 
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pass through the bubble wall they (1) evolve in 
space and time according to our equations of motion 
and do so coherently with the Higgs field to which 
they couple; and, (2) populate the bubble preserving 
the number density to the extent allowed by overall 
energy-momentum conservation, as discussed in the 
appendix. 

Secondly, it is clear that the magnetic field ex- 
tends more deeply into the collision region in our 
theory, which makes the volume-averaged azimuthal 
field in the surface of the ring even larger than it is 
in the Abelian Higgs model. This possible source of 
enhancement is however mitigated by the fact that 
the finite conductivity a damps the magnetic field 
in the interior of the ring. The damping is effective 
beginning at depth Sp r i ng from the surface of the 
ring, 



where 



$ Pring ~ 2 



6.7T C 



(61) 



(62) 



The ring expands with the time elapsed after the 
initial collision 5t as 



Pr 



VmuSt , 



(63) 



where R is the radius of a bubble when the collision 
occurred and v is the wall speed [Tl| . 

Note also that the magnetic field in our theory 
drops to zero at the outer boundary of the expanding 
overlap region whereas it terminates suddenly in the 
Abelian Higgs model. This unphysical feature was 
discussed in Ref. [l3[ and was traced to the abrupt 
change in the boundary condition at the point of 
collision. The same behavior occurs in our theory in 
the fields, but its effect on the magnetic field 
is evidently mitigated when we evaluate the current 
with the charged as its source. 

Next, we want to examine the magnetic field as 
we move along the z direction at fixed radial dis- 
tance p from the collision axis. By comparing the 
results shown in Fig. [5] and Fig. [5] for our theory and 
that of the Abelian Higgs model, respectively, we 
again see similar magnitudes of the magnetic fields. 
The azimuthal magnetic field in the Abelian Higgs 
model again reaches its peak at the surface of the in- 
tersection region, whereas in our theory is generally 
small there. Thus, again we see that when we take 
the source of the electromagnetic current to be the 
charged fields, the magnetic field created in col- 
lisions behaves in a smoother fashion than it does in 
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FIG. 5: Magnitude of the azimuthal magnetic field cal- 
culated in the theory of the present paper. The field is 
shown as a function of distance z along the axis of col- 
lision at a distance p = 1 from the axis of collision for 
times St = 5, 10, 15, and 20. 



0.04- 

0.03- 
0.02- 
0.01- 



W 1 



FIG. 6: Magnitude of the azimuthal magnetic field cal- 
culated in the Abelian Higgs model. The field is shown 
as a function of distance z along the axis of collision 
at a distance p — 1 from the axis of collision for times 
St = 5, 10, 15, and 20. 



the Abelian Higgs model with step-function bound- 
ary conditions. The tendency to peak at the edge 
of the intersection region makes strong longitudinal 
oscillations in B in the Abelian Higgs model, partic- 
ularly when we look deep within the overlap region. 
Because this region characterized by relatively small 
p has had more time to evolve than the regions closer 
to the outer ring where p is large, the conductivity 
has its most pronounced effect, strongly damping 
the fields as well as concentrating them into a single 
peak [12j . 

The quantity we need to obtain from our calcu- 
lation for determining the effectiveness of magnetic 
field creation during the electroweak phase transi- 
tion seeding the galactic and extragalactic magnetic 
fields we see today is the azimuthal magnetic field 
strength B(R) created at the surface of the ring that 
is formed by the colliding bubbles. In the Abelian 
Higgs model, the value of B{R) is given by Eq. (24) 



13 



of Ref . [l2j as a function of v wa u , the radius of the 
bubbles R at the time of collision t = t C) and the 
radius p r i ng at the time of completion of the phase 
transition. This result does not depend explicitly on 
the conductivity a. This is because the conductivity 
damps out the magnetic field in the interior of the 
bubble but not at the outer surface of the ring where 
it is the largest. However, the finite wall speed does 
damp the magnetic field at the surface as shown in 
Eq. (24) of Ref. fl2| . Because this damping is a kine- 
matic effect, it scales with the value of the magnetic 
field at the surface of the ring. Hence, from knowl- 
edge of the value of the magnetic field in the ring, we 
can obtain the appropriate B(R) for our theory by 
scaling it relative to the value of the magnetic field 
in the ring in the Abelian Higgs model. It is clear 
from Fig. [3 Fig. II Fig. El and Fig. [U that this scal- 
ing indicates that our seed fields will lead to about 
the same magnetic field as the estimate in Eq. (30) 
of Ref. [l2j for the galactic dynamo. 



VII. SUMMARY AND CONCLUSIONS 

Methods suitable for exploring magnetic seed field 
creation during a first-order EW phase transition 
have been developed for the MSSM. We obtain equa- 
tions of motion from the EW Lagrangian, omitting 
fcrmions. One of the most important consequences 
of these equations is that the fields of the non-abelian 
sector are the only charged ones that appear in the 
electromagnetic current, and therefore these fields 
(the charged gauge bosons W ± ) are the sole source 
of the magnetic seed fields. This is in sharp contrast 
to the Abelian Higgs model ([ll|, [13] , [3), where 
the magnetic field arises from gradients in the phase 
of the Higgs field. 

We studied numerically the magnetic fields pro- 
duced in collisions of the bubbles of the first-order 
electroweak phase transition with our equations of 
motion and found, for gentle collisions (where fluc- 
tuations in the Higgs field are relatively small in 
bubble collisions), that the W fields required for 
calculating the current may be obtained by solving 
linear equations. Numerical solutions of the par- 
tial differential equations were obtained by applying 
simple boundary conditions that fix the fields at the 
time of bubble collision, similar to those employed 
in Refs. [Ill, d3, and [Hj]. The magnetic seed fields 
we obtained with the charged as their source are 
comparable to those estimated in the Abelian Higgs 
model. 

These results establish the electroweak phase tran- 



sition in the MSSM as a promising source for pro- 
duction of seed fields for large-scale galactic and 
extra-galactic fields observed today. In our future 
work [211 ] on magnetic field generation we will ex- 
plore sensitivity to some of the more restrictive as- 
sumptions we make in this work. 
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Appendix: Normalization of Charged 
Fields 

The process by which thermal modes of the 
plasma make the transition into the bubble is a com- 
plicated many-body problem that we are unable to 
simulate. In order to obtain numerical estimates, we 
introduce instead a model for determining the nor- 
malization of the charged VF ± fields. Our main as- 
sumption is that the number density of W quanta 
in bubble is the same as the number density of the 
quanta in the thermal plasma in the modes from 
which they originated. The density of the W out- 
side the bubble that are able to make the transition, 
pw, is 



Pw(T) = 



^ d 3 k 



(2tt) 3 e k / T 



(64) 



where k is fixed by overall energy-momentum con- 
servation. Energy conservation is an issue because 
in the plasma the are massless (having energy 
hk), whereas inside the bubble their energy is the 
same as their mass (see Eq. (|36|) ). Overall energy- 
momentum conservation is a weaker constraint than 
insisting that energy and momentum be conserved in 
each collision, and is a valid constraint when multi- 
ple simultaneous collisions occur at the bubble wall, 
as we assume to happen at the high densities of inter- 
est here. Application collision-by-collision will lead 
to a finite wall speed |2(|. The total energy of the 
in a single bubble is given by their mass times 
the density in Eq. (fS4"]l . 

Overall energy conservation thus determines the 
portion of the spectrum (fc < k) that populates the 



14 




100 120 \140 160 180 GeV 



FIG. 7: Spectral limit 



bubble as it grows, 
mw 2 

= 2 



h ^ d 3 k 1 



( 27r )3 e k/T _ I 

d 3 k i 



k 



(65) 



o (2tt) 3 e k / T - 1 



The factors of 2 arise from the fact that there are 
only two spin directions for massless vector bosons 
in the plasma. We show the solution k(T) of Eq. ([65]) 
in Fig. [7J In the region of temperature shown, k 
3rnn//2. The density pw of W ± inside the bubble 
is then 



Pw{T) = — \ 
71-2 Jo 



T 3 f k ^' T x 2 dx 



(66) 



The number of quanta filling the bubble clearly 
grows as the bubble expands and displaces a larger 
volume of the plasma. Thus, inside the bubble the 
number of W boson pairs N is N = pw{T)Q, 
where fl = 47r/3(i?o + i 2 ) 3/2 is the volume of the 
bubble. 

For a single , assuming relativistic normaliza- 
tion of wave functions, 



w = 



(67) 



where is the volume of the bubble. Therefore, for 
N W bosons in a bubble, 



w = VN/\/2m w ri 



' pw(T) 
2mw 



(68) 



For k given in Fig. [Jj w is nearly linear in temper- 
ature and is given over the same temperature range 
as 



W (T) 2 = -127. + 4.49T 



(69) 



in units of GeV. It is unlikely that more W from 
the plasma than allowed by energy-momentum con- 
servation would drop into the coherent mode since 
the excess energy would lead to heating that would 
impede the phase transition. 
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